On Size and Shape of the Average Meson Fields 
in the Semibosonized Nambu &; Jona-Lasinio Model 



R. Wunsch\ K. Goeke^, Th. Meissner^ 



^ Institut fiir Kern- und Hadronenphysik, Forschungszentrum Rossendorf e.V. 
Postfach 51 01 19, 01 314 Dresden, Germany 

^ Institut fiir Theoretische Physik II, Ruhr-Universitat Bochum, 
Postfach 10 2148, 44 780 Bochum, Germany 



3 



Institute for Nuclear Theory, University of Washington, IIN-12 Seattle, 
WA 98195, USA 



Abstract: We consider a two-flavor Nambu & Jona-Lasinio model in Hartree ap- 
proximation involving scalar-isoscalar and pseudoscalar-isovector quark-quark interac- 
tions. Average meson fields are defined by minimizing the effective Euklidean action. 
The fermionic part of the action, which contains the full Dirac sea, is regularized within 
Schwinger's proper-time scheme. The meson fields are restricted to the chiral circle and 
to hedgehog configurations. The only parameter of the model is the constituent quark 
mass M which simultaneously controls the regularization. 

We evaluate meson and quark fields self-consistently in dependence on the constituent 
quark mass. It is shown that the self-consistent fields do practically not depend on the 
constituent quark mass. This allows us to define a properly parameterized reference field 
which for physically relevant constituent masses can be used as a good approximation to 
the exactly calculated one. The reference field is chosen to have correct behaviour for 
small and large radii. 

To test the agreement between self-consistent and reference fields we calculate several 
observables like nucleon energy, mean square radius, axial-vector constant and delta- 
nucleon mass splitting in dependence on the constituent quark mass. The agreement is 
found to be very well. 



1 Introduction 



The model of Nambu & Jona-Lasinio (NJL) has been used quite successfully as effective 
chiral theory for low and medium energy hadronic phenomena. First it has been applied to 
vacuum and meson properties as well as medium effects (for reviews c. f. [§, ^ Q). Later 
on it turned out that also baryonic systems (nucleons and hyperons) can be described 
within this model (for a review c. f. f^). Starting from a semi-bosonized version |P, |^ 
with scalar-isoscalar and pseudoscalar-isovector interaction and treating the meson fields 
classically various authors have shown that for constituent quark masses M ~ 350 MeV 
it is possible to get self-consistent solitonic solutions with baryon number B=l consisting 
of 3 valence quarks in addition to the polarized Dirac sea P]-|T^. Because of the non- 



renormalizability of the Nambu & Jona-Lasinio model the sea contribution diverges and 
has to be regularized. The parameters of the model can be fixed to the physics of the 
meson and vacuum sector, mainly to the weak pion decay constant and the pion mass 
m,r- In doing so only one free parameter remains open for the baryonic sector for which 
we take the constituent quark mass M. 

The self-consistent determination of the meson fields is a time-consuming numerical 
procedure. Changing the parameters of the model or the regularization scheme the proce- 
dure has to be repeated again. So it might be helpful to look for an analytic parametriza- 
tion of the selfconsistent profile function G(r) of the solitonic solution which approximates 



the exact 0(r) as well as possible. Within the restrictions to hedgehog configurations fTS 



and to the chiral circle the meson fields are uniquely described by the profile function 
0(r). In the course of our calculations we noticed to a very large extent an independence 
of this profile function on the constituent quark mass M. It is the aim of this paper to 
investigate this dependence quantitatively and to look for a general function which may 
approximate the profile function, if possible independently of M. 

In section 2 we review the main ideas of the semi-bosonized and regularized Nambu & 
Jona-Lasinio model for two flavors and introduce observables characterizing the quark and 
meson configuration. The dependendence of the self-consistent meson profiles 0(r) on the 
constituent quark mass M is investigated in sect. 3 within a wide range (350 MeV < M < 
1000 MeV). Here we compare the self-consistent profile with a reference profile Q^'^^{R; r) 
given by a simple arithmetic expression obtained from an asymptotic expansion of the 
equation of motion at r ^ and r — >■ cxd. In sect. 4 we evaluate several observables like 
nucleon mass, mean-square radius, axial-vector coupling constant and delta-nucleon mass 
splitting using both the self-consistently determined profiles and the standard profiles. 
The comparison of both values tests and illustrates the quality of the reference profile. 



The regularized and bosonized Nambu & Jona- 
Lasinio model and its observables 



The details of the following section can be found in ref. ||^, 16, 17, 21|. Here we shortly 
review these parts of the formalism which make this paper self-contained. 
We consider a two-flavor NJL lagrangian 



NJL 
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for the quark fields q{x) {u and d quarks of = 3 colours). Here G is the coupling 
constant, r the vector of Pauli-matrices and m the average current mass of the light 



quarks m = (m„ + m(i)/2. The theory with only quark degrees of freedom is converted 
into an effective quark-meson theory by means of standard path-integral bosonization 
0, |T9|. The mesonic fields are restricted to hedgehog configurations and to the ciral 
circle. The resulting semi-bosonized theory is described by an effective (Euklidean) action, 
whose dynamical quantity is the profile function 6(r) [|10], |1^ 
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consists of a quark part 



The Aeff 
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with the quark hamiltonian 
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and of a meson part 
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where T is the total Euklidean time interval. The quark-meson coupling strength g 
in the hamiltonian (^) and the parameter A in the mesonic action @ are related to 
the interaction strength G in the NJL lagrangian (|I]) via G = g'^/X'^- Since we are 
interested in stationary field configurations we have assumed the meson fields to be time- 
independent. Hence the eigenvalues Ea of the hamiltonian (^) are real. The vaccum state 
(with broken chiral symmetry) has 0(r) = and and the eigenvalues and eigenf unctions 
of corresponding hamiltonian are denoted by and $^(r), respectively. The symbol Sp 
indicates functional and matrix (spin, isospin, colour) trace 



SpO = Nctr-ytr-r J d'^XE {xe\0\xe) 
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with the Euklidean space-time vector = (t, r) and its volume element d^XE- The 
effective quark action (|]) can be divided into a valence contribution 



Aaim = A[0](M)-^n0](-" = o) 

and contribution 
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The valence contribution (|^) depends on the chemical potential fi which will be adjusted 
such that the resulting soliton has baryon number one. This can be achieved by a chemical 
potential fi > which is slightly larger than the energy e^ai of the lowest quark level with 
positive energy (valence level). However, the meson field may be so strong that the valence 
level is bent down to the Dirac sea [e^ai < 0). In this case the Dirac sea has one level 
more than in vacuum state and therefore carries baryon number one. Describing a baryon 
with baryon number one we have to choose = 0, in this case, making the Fermi sea 
(0 < < /i) empty. 



The sea contribution diverges and is regularized. Applying Schwinger's proper- 

, where A is the regularization 
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parameter, we get a regularized sea contribution 



Here the regularization function is given by 
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with the incomplete Gammafunction Y[x^a). The total regularized effective action reads 



The profile function G(r) is determined by 

[6] 
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yielding the equation of motion 
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The valence contributions S^aiif^) and Pvaii'f') are obtained for the valence level a = val. 
They vanish if the valence level dips into the Dirac sea {e^ai < 0). The $0,(7") are the 
eigenfunctions of h. The sea contributions have been regularized within the proper-time 
scheme with a regularization function 



i?„(£,A) = ^ / dtr^'^e-'' 



sign{e) /l 

I2' A2 



1/A2 



TT 



The parameters of the model will be fixed by adjusting the vaccuum and meson sectors 
as described in detail in ref. [0]. The only remaining free parameter is the constituent 
quark mass 



M = gU 

which will be varied within reasonable hmits. 

In this paper we consider the (isoscalar) baryon density 
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The colon in the matix element of eq. (|20|) indicates subtraction of the corresponding 
value calculated in the vacuum state. Like for any other observable the valence contri- 
bution (0) vanishes if the Fermi sea is empty. The Dirac sea contribution (^) has not 
been regularized because it is finite and fulfills the exact normalization condition (baryon 



number zero/one) only in the unregularized case |jT6[. The size of the density distribution 
is characterized by the mean square baryon radius 
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Another quantity characterizing a quark configuration is the axial density 



It consists of a valence contribution 
and of contribution 
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which has been regularized with the regularization functions i?m(^; A) defined in eq. (|T8D 



21| . The axial density determines the axial- vector coupling constant of the proton 

-2 /rfVAJr), 



gA = --i ia"rAo[r), (27) 
where an additional factor (—1/3) is incorporated which results from the projection onto 



the isospin quantum number T = 1/2 of the proton |23 



The total energy E' of a static quark-meson configuration is given by 

E = El, + Et^'3 + ^ ^csB 
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with the valence-quark energy 



rpq 



1 

r 



d 



and regularized sea-quark contributions 
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The meson energy is given by 
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For time-independent fields the classical meson and sea-quark energies differ from the 
corresponding effective actions P) only by a factor T. An additional contribution to 
the energy results from /i dependence of the valence-quark action (|^. 

Finally we consider two spurious contributions to the quark energy which result from 
the mean- field approximation and have to be subtracted from the total energy (pSl) . As 
shown in p, 22 1 the static hedgehog contains a center-of-mass motion with the energy 



: Jd-^rq'{r) (-v ) q{r) 



EcMM - - ^ , (32) 

where E is the rest mass represented by the total hedgehog energy (|^). The expectation 
value of the square of the total quark momentum P consists of a valence- and a sea- 
contribution given by 



val 
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Another correction term results from the quark rotational degrees of freedom |23, 24, 25 



It is described by a moment of inertia J, which can be calculated within the semiclassical 
cranking approach p6| and consists of a valence contribution 
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and of a regularized sea contribution 
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Within the proper-time scheme the regularization function is given by 
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In the limit A ^ oo one gets the well-known Inglis formula |^ for the moment of inertia. 
The incomplete error-function is given by er/c(x) = dt. Since the energy 

correction for the nucleon and the A isobar are different the moment of inertia gives rise 
to a mass splitting between both particles which is given by p5 
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3 The self-consistent meson profile and its depen- 
dence on the constituent quark mass 

The equation of motion (|1^) is an implicite and nonlocal equation for the profile function 
0(r). Because of the dependence of the expectation values S{r) and P(r) on the eigenval- 
ues and eigenf unctions of the hamiltonian h the right-hand side of the equation of motion 
is a functional of the profile function. We determine B(r) iteratively. Starting from a 
reasonable profile 0°(r) we determine eigenfunctions $q(t*) and eigenvalues £° by diag- 
onalizing the hamiltonian (^) within an appropriate basis. By means of the equation of 
motion (|13D and the auxileary functions (|l4| - p!^ ) we get an improved profile function B^(r). 
Continuing this procedure to convergency we get the self-consistent profile function. 

We represent the eigenfunctions of the hamiltonian within a discrete basis introduced 
in [^. It is defined within a spherical box with a radius D which is several times larger 
than the extension of 6 field. Details of the procedure are described in [|1^, 0, [1^ . 

For all our calculation a box radius D = 15 /M has turned out to be sufficiently large. 
The discrete basis was limited by a maximal momentum K^ax = 8 M, which is more than 

4 times larger than the corresponding regularization parameter A. 

We have numerically determined the self-consistent profiles for 13 values of the con- 
stituent quark mass M between 350 and 1000 MeV. As illustrated in fig. 1 the profiles are 
nearly independent of the mass parameter. None of the calculated profiles leaves the nar- 
row corridor limited be the two broken lines. At masses lower than 350 MeV the iteration 
converges to the vaccum field (G(r) = 0). Here the corresponding interaction strength G 
is not strong enough to keep a solitonic configuration together. 

The self-consistent profiles can be approximated by a reference profile 



with the matching point 
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Rm = (l + + ^li^nR?] ) (40) 



and an empirically determined radius parameter 



R = 0A2fm. (41) 



The reference profile (|39|) with the matching point (^Of ) interpolates smoothly between 
the correct asymptotic behaviour at r — and r — > oo following from the asymptotic 
expansion of the equation of motion (|13|) . It constitutes - at least visually - a fair approxi- 
mation of the self-consistent profiles obtained after a time-consuming iteration procedure. 
Further tests of the reference profile will be performed in the next section. 



4 Testing the reference profile on nucleon observab- 
les 

In sect. 2 we considered several expectation values characterizing a quark or meson config- 
uration. All these quantities are functionals of the profile function within our model. We 
have tested the quality of the reference profile (|39|) and evaluated the observables using 
both the reference and the self-consistently determined profiles. Fig. 2 shows the total 
energy (|28| ) and its components (p9| , |30| , ^1]) , the mean square radii (^3]) of the baryon 
density, including their valence- and sea-quark contributions, and the axial-vector cou- 
pling constant of the proton (pTj) calculated with either profile. The kinks in the valence 
and sea contributions at the critical constituent mass M = M^ip ~ 750 MeV result from 
the valence level which leaves the Fermi sea and joins the Dirac sea. The behaviour of 
the regularization functions (0, |18], at e — guarantees that the sum of valence 



and regularized sea contributions is a smooth function of the constituent quark mass M. 
Fig. 2 illustrates that nicely. 

The only noticeable difference between the values for self-consistent and reference 
profiles appears in the valence and sea contributions in the vicinity of MfUp. For the 
self-consistently determined profiles, the valence level dips into the Dirac sea aX M 
IbQMeV. This point is shifted to M 725 for the reference profile. The deviation 
is another evidence for the more sensitive dependence of valence and sea contributions 
on details of the profile function, while their sum is quite insensitive. One should note, 
however, that the physically relevant region for the constituent mass is around M = 
400 MeV, where the nucleon observables get reproduced by the reference profile very 
well. 

The calculated nucleon observables are in sufficient agreement with similar calculations 



16, 21, 17]. The too small value of the axial-vector coupling constant {gA ~ 0.6 ~ 0.8) 
in comparison to the experimental value {q'a^ ~ 1.25) is a lack shared by many chiral 
models of the nucleon. However it is rather the aim of this paper to compare between 
two theoretical approaches than to reproduce the experimental values. 

To complete our test of the reference profile we evaluate corrections to the static 
hedgehog configurations due to zero-point modes and spin-isospin quantization. Fig. 3 
shows the energy Eqmm (0) of the center-of-mass motion, the moment of inertia / 
including the components ( ^5]) and (^61), and the resulting nucleon-delta mass splitting 
We establish an excellent agreement between the values for both kinds of profiles 



in the physically relevant mass region below 600 MeV. Larger deviations appear at M ~ 



600 MeV. Here the energy corrections are so large that the perturbative approach used 
for their determination is already not justified. 



5 Conclusions 

We have self-consistently calculated average meson fields for the SU(2) Nambu & Jona- 
Lasinio model with scalar-isoscalar and pseudoscalar-isovector couplings in Hartree ap- 
proximation. The fields are restricted to the chiral circle and to hedgehog configurations. 
Infinite quark contributions are regularized within Schwinger's poper-time scheme. 

The numerically determined meson profile functions turn out to be nearly independent 
of the constituent quark mass. They can be approximated quite well by a unique reference 
profile given by a simple arithmetic expression, which interpolates between the correct 
asymptotic behaviour at large and small radii. It is shown that the reference profile 
does not only approximate the self-consistent profiles but also reproduces the relevant 
observables of the quark and meson configurations. 

We conclude that many of the properties of the Nambu & Jona-Lasinio lagrangian 
can be studied using the reference profile instead of applying the time-consuming deter- 
mination of the self-consistent profile. Changing the constituent quark mass M only the 
strength g of the quark-meson coupling is infiuenced {g ~ M), while shape of the meson 
fields is almost independent of M. The absolute strength of the meson fields was fixed 
by the restriction to the chiral circle which can be justified from an extended NJL model 
implementing the trace anomaly of QCD ^ . If an accurate determination of the 



self-consistent profile is necessary, the reference profile may serve as a suitable starting 
profile. 

The reference profile can be compared with the Woods-Saxon potential describing the 
average field inside an atomic nucleus. Most of the nuclear properties are sufficiently well 
described by this model potential which is rather determined by a simple ansatz (with 3 
parameters) than by a sophisticated Hartree or Hartree-Fock procedure. 
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FIGURE CAPTIONS 



1. Self-consistet profiles in the mass region 350 MeV < M < 1000 Me V" All self- 
consistently calculated profiles fit in the area marked by the broken lines. 

The full line represents the reference profile ( ^9| ) with R = 0.42 fm. 

2. Nucleon observables in dependence on the constituent quark mass M calculated 
with self-consistently determined profiles {full lines) in comparison to the reference 
profile B^^'^(-R;r) defined in eq. (p9|) with R = 0.42 /m. {broken lines) 

Upper part: Total energy E (|28| ) and its components E^ai (|29D, E^J'J^ (|30| ) and E"^ 

m- 

Central part: mean square radius R ( ]23|) and its valence {Rvai) and sea contributions 
{Rsea) calculated with the corresponding densities (^Tf) and respectively. 
Lower part: Proton axial- vector coupling constant qa (total value only). 

3. Energy corrections to the static hedgehog energy calculated with self-consistently 
determined profiles {full lines) and with the reference profile {broken lines) in de- 
pendence on the constituent quark mass M. 

Upper part: Energy of the center-of-mass motion Ecmm dH) . 

Central part: Moment of inertia I as the sum of the valence part I^ai (0) and the 

sea contributions I sea (|36|) . 

Lower part: Delta- nucleon mass splitting (|38|). 
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